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ABSTRACT 

We  present  here  an  0(0^)  deterministic  and  an  0(n)  probabilistic 
algorithm  for  computing  the  area  of  the  union  of  n  given  circles.  The 
accuracy  of  the  probabilistic  algorithm  is  controlled  by  the  algorithm 
iraplementer  (it  can  be  set  up  to  any  desirable  value).  The 
probabilistic  algorithm  is  optimal  in  the  sense  that  0(n)  time  is 
needed  to  compute  the  sum  of  n  real  numbers,  under  quite  general  models 
of  computation.  The  probabilistic  linear  algorithm  can  be  extended,  in 
a  straightforward  manner,  to  compute  the  volume  of  the  union  of  many 
similar  objects  (e.g.  spheres)  in  k  dimensions.  Its  time  complexity 
is  then  0(nk)  . 

Our  proof  of  fast  convergence  of  the  probabilistic  method 
benefited  a  lot  from  a  recent  work  of  R.  Karp  and  M.  Luby  on 
reliability.  Our  technique  was  also  motivated  by  Karp  and  Luby's  work, 
and  some  of  our  notation  is  heavily  affected  by  their  notation. 
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1 .   Introduction 

The  problem  of  estimating  the  area  of  the  union  of  many  circles  in 
the  plane  was  first  posed  by  [Sharaos,  78].  We  present  here  an  0(n  ) 
deterministic  and  an  0(n)  probabilistic  algorithm  for  solving  the 
problem.  The  recent  work  of  [Karp,  Luby ,  83]  on  estimating  the  failure 
probability  of  an  n-component  system  helped  us  a  lot  in  stating  the 
algorithm  and  formulating  the  proof  of  fast  convergence.  The 
probabilistic  algorithm  can  be  modified  to  compute  the  area  of  the 
union  of  other  planar  objects  (of  fixed  description)  e.g.  triangles  or 
boxes.  The  time  is  again  0(n).  The  algorithm  can  also  be  extended  to 
compute  the  volume  of  the  union  of  n  spheres  in  k  dimensions,  in  time 
0(nk).  The  algorithm  has  similar  time  complexity  for  k-dimensional 
objects  other  than  spheres,  provided  that  each  object  has  a  fixed 
description  (not  dependent  on  n) . 

A  very  recent  work  of  [Sharir,  83]  can  be  extended  to  provide  an 
0(n  log  n)  deterministic  algorithm  for  computing  the  area  of  the  union 
of  n  circles  in  the  plane.  The  idea  does  not  generalize  easily  to  3 
dimensions.  No  efficient  algorithms  for  the  problem  of  computing 
volumes  of  unions  of  objects  in  more  than  2  dimensions  were  presented 
in  the  past . 


2.    The  problem  in  2  dimensions 
2.1   General  remarks 


The  input  is  a  list  of  n  triples  (x^^,  y.,  r.),  I  =  l,...,n,   where 
(x^,   y^)  are  the  Cartesian  coordinates  of  the  center  of  the  circle  C, 
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n 
of  radius  r^  .   The  problem  is  to  compute  the  area  of  the   union  U   C,  . 

Let   U   be   the   area   of   the  union.   Let  E(C.)  denote  the  area  of  the 

n 

circle  C .  .   Let  ),  =  ),   E(Cj^).   The  part  of  the   plane   covered   by   the 

i  =  l 
union  of   the   n  circles  is  partitioned  into  nonoverlapping  segments, 

whose  boundaries  are  circular  arcs  (pieces  of  the  circumferences  of  the 

circles   C.).    Let  s^,...,Sj^  be  an  enumeration  of  these  segments.   Let 

ra 

E(s.)  represent  the  area  of  segment  Sj^.   Clearly,  ^   =    \     E{s^).      In  our 

i=l 
algorithms,    we   assume   that   all   calculations   are   done   in  exact 

arithmetic . 


Remark  2.1    The  number  of  segments,  ra,  is  O(n^) 

Proof:  Let  us  assume  we  have  already  i  circles  drawn  and  we   next   draw 

the   i+1    circle.   In  the  worst  case,  2i  new  points  of  intersection  of 

circumferences  are  introduced  (2  for  each  of  the  i   circles).    Let   us 

call   them  vertices.   So,  the  total  number  of  vertices  for  n  circles  is 
n 

v  <  )_   2i  =  0(n  ).   If  we  take  out  degeneracies,  a  vertex  is  shared   by 

i  =  l 
at   most    4    segments,    hence   the   total   number   of   segments   is 

m  <  4v  =  0(n2)  .  • 

The  above  remark  leads  to  a  straightforward  0(n  )  time  and  0(n  ) 
space  algorithm  for  computing  U.  This  algorithm  keeps  a  list  of  the 
segments  created  by  the  i  first  circles.  Each  segment  is  kept  as  a 
sequence  of  circular  arcs,  this  sequence  being  the  arcs  which  surround 
the  segment,  listed  in  the  counterclockwise  direction- 
Each   circular   arc   is   a   triple  of  a  start  point,  an  end  point  and  a 
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curvature.   Clearly,  a  segment  created  by  i  circles  can  have  at  most   i 
surrounding  arcs. 


2.2   An  0(n2)  algorithm 

This  algorithm  keeps,  in  each  stage,  just  a  list  of  circular  arcs, 

surrounding  the  current  union  (of  the  first  1  cycles).   The   algorithm 

is  as  follows: 

Algorithm  1 

Initially:   i  =  1.    Let  iil)      be   the   circumference   of   C, .   Do  the 

following  stage  for  i  =  2  to  n. 
/*  Begin  stage  i 

Let  Z(i-l)  be  the  current  list  of  arcs  of  the  union  (The  arcs  of  each 
connected  component  are  listed  in  counterclockwise  direction).   Find 
its  intersections  with  the  i*-"  circle  (There  can  be  0(i)  such 
intersections).   Updatae  the  list  by  removing  arcs  which  are  now 
interior  to  the  new  union  and  by  adding  the  surroundings  of  the 
part  of  the  i+1    circle  which  is  not  covered  by  the  union  of 
the  i-1  previous  circles.   The  new  list  is  the  list  ^(i). 

/*  End  stage  i. 

Finally,  compute  the  area  enclosed  into  the  arcs  of  the  list  5,(n). 

Lemma  2.1  Algorithm  1  takes  0(n2)  time  and  0(n)  space. 

Proof :  The  Algorithm  needs  0(i)  time,  to   update   the   list   2.(1-1)   at 
state   i,   since   £(i-l)   may  have  up  to  L-1  arcs,  one  per  each  circle. 
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So,  the  total  time  spent  during  the  n  stages  is  0(n  ).  The  area 
calculation  takes  0(n)  time.  (The  area  calculation  can  be  done  by 
computing  0(n)  integrals,  one  for  each  arc,  and  adding  them  with  a  plus 
or  minus  sign).  The  space  needed  is  that  of  storing  the  list  of  arcs 
which  surround  the  union.  This  list  cannot  have  more  than  n  arcs.  • 
Note.  Algorithm  1  computes  not  only  the  area  of  the  union  but  also  the 
boundary  of  the  union  of  the  n  circles. 


2.3    A  straightforward  Monte-Carlo  approach. 
2.3.1   Definitions  -  Terminology 

Following  [Karp,  Luby,  83],  let  U  be  the  estimate  of  the   area   of 
the  union  (U) ,  produced  by  a  randomized  (Monte  Carlo)  algorithm  MC . 

Definition.  MC  is  called  an  (e,6)  algorithm  of  time  f(n),  if 
Ve,6  €(0,1),  one  can  select  a  running  time  f(n)  (n  =  number  of  circles) 
of  MC ,  such  that  the  estimate  U  produced  has  the  property 


ProbilLi^l    >   e|    <    5 


(See  [Karp,  Luby  83]  for  a  similar  definition)  If  f(n)  =  0(n),  MC  is  a 
linear  (e,5)  algorithm. 

Definition.  Given  the  set  of  n  circles  C.  ,  i  =  1  ,  ,  .  . ,n  In  the  plane, 
and  a  point  P,  let  cover(P)  be  the  number  of  circles  to  which  P  belongs 
(we  consider  circles  to  be  closed  sets). 
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2.3.2.   The  straightforward  approach 

Let  us  assume  that  we  can  find  (in  time  a(n))   a   closed   curve   A 

(enclosing  area  E(A)),  which  completely  encloses  the  n  circles.   Assume 

also  that  we  have  a  fast  way  (of  time   T(n)   per   point   selection)   of 

selecting   points  within  A  (or  the  interior  of  A)  in  a  uniformly  random 

manner  (meaning  that  the  probability  of  a  point  selected  to   fall   Into 

an  elementary  area,  Ae  ,  around  a  specific  point  Pq ,  is  the  same  VP^   (A 

plus  interior  (A))  and  it  is  equal  to  •  Ae.) 

E(A) 

Finally,  let  3(n)  be  the  time  it  takes  to  decide  if  a  given 
(arbitrarily  selected)  point  P  belongs  to  the  union  of  the  n  circles. 

Straightforward  MC 

1     E(  A) 

(I)  We   select   N  >  [ -   ij  points,  uniformly  randomly. 

(Note  that,  to  do  this  selection  we  don't  need  to  know  U,  just  an  upper 
bound  on  E(A)/U  is  sufficient). 

(II)  For  each  point  selected  in  (I),  we  test  if  it  belongs  to  the 
union  of  circles.  Let  M  <  N  be  the  number  of  points  found  to  belong  In 
the  union. 

(III)  Output  U  =  ^  .  E(A)  . 

N 

Lemma  2.2.   The  straightforward  MC  method  is  an  (e,5)  algorithm  of  time 
0(a(n)  +  T(n)B(n)  (i^  -  1)  ^J  . 


Proof .   Let  x.  =  1  iff  he  i*"   selected  point  belongs  to  the  union,   and 

N 

x^  =  0  else,  for  i  =  l,...,n.   Then  Xj+...+Xf^  =  M  and  mean(M)  =  \      mean 

i=l 
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(x.)  =  N  .  -^.      So,  mean   (U)  =  ^i^  niean(M)  =  ii^  N  J^  =  U(*); 
^        E(A)  N  N       E(A) 

hence  U   is   unbiased-    Let  a        =   variance   (x^)   i   =  1,...,N.   By 

?     U        U 

elementary   calculation,   o^  =       -  (**)•    Then,    variance 

E(A)     £2(^) 

(M)  =  a   •  N  hence 


variance(U)  =  ^  E^CA)  =  ^^Il1!1  =   ^(AIlHl^ 
n2  n        N 


by  (*),  (**). 

By  Chebyshev  Inequality  then, 


ProbilU-U|  >  eU}  <  variance(U)  ^  E(A)U-U   _  _1_  ^E(A)_  _  ^^ 

ir2ii2  vTi:.2n2  Mc-2  U 


Hence   Probj-! >   e}    <    6    ,  QED 


Clearly  the  running  time  is  a(n)  (to  select  A)  plus  0(N«T(n) • S(n) ) .   • 

2.4   How  to  improve  on  the  straight orward  MC . 
(An  O(nlogn)  (£,5)  algorithm). 
A  reasonable  value  of  a(n)  is  0(n)  (Eg.   we  can  arbitrarily  select 
a  center  and  then  find  the  most  distant  of  the  centers  of  the   C.'s   to 
that   center,   add   to  the  distance  found  the  maximum  radius  and  draw  a 
circle  for  A) . 

If  we  don't  allow  preprocessing,  then  3(n)  =  0(n),   leading   to   a 
time  complexity  of  0[n«T(n)  •  (:i^^ — L-    1)J.   It  should  be  clear  that  if 

(1)  A  is  easily  describable  so  that  T(n)  =  0(1)  and 

'C'  /    K   \ 

(2)  — _ —  =  0(1),  then  we  would  get  a  linear  (£,6)  algorithm. 
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Requirement   (2)  implies  that  E(A)  is  close  to  U.  One  such  A  would 

n 
be  the  boundary  set  of  the  union  U   C.  itself,  but   then   the   task  of 

1  =  1   ^ 

selecting  points  from  the  interior  of  A  uniformly  randomly  becomes  a 
hard  task,  since  a  point  belonging  to  k  >  1  circles  cannot  be  counted 
as  a  "whole  point"  when  computing  M. 

A  first  modification  to  MC  (Call  it  Extended  MC ,  EMC)  is  the 
following: 

Extended  MC 

(1)  Let  A  be  the  boundary  of  the  union  of  the  n  circles.  (i.e 
E(A)  =  U) 

(2)  To  select  a  random  point  in  the  interior  of  A  we  select  one  of 
the  C.  at  random  (Call  it  C. )  and  then  we  select  one  random  point  P^  of 

the  circle  selected. 

N 

(3)  We  compute  M  =  )^   (cover  (P.))~ 

j  =  l  ^ 

M 

(4)  We  output  U  =  ^ 

N 

It   is   relatively   easy   to   prove  that  if  N  >  then  EMC  is  an 

(e,5)  algorithm.  Its  complexity  depends  on  the  time  F(n)  to  find  cover 
(?)  for  each  P.  By  using  the  Voronoi  diagram  construction  as  a 
preprocessing  idea,  we  can  compute  cover  (P)  in  time  O(log  n)  for  a 
given  P.  This  will  result  in  an  0(n  log  n)(e,5)  EMC  algorithm. 

In  the  next  section  we  will  improve  upon  the  time  complexity  of 
the  0(n  log  n)  EMC  by  using  the  following  two  probabilistic  ideas: 

(a)  We  shall  assign  greater  probability  of  selection  to  points 
belonging  to  circles  of  bigger  area. 
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(b)   We  will   estimate   covers   through  a  very  fast  probabilistic 
experiment . 


3.   A  linear  (e,5)  algorithm. 

3.1   Linear  Monte  Carlo  (LMC)  Algorithm  Initialization: 

c  is  a  small  constant  >  20 

S:=0  ,  number  .of .trials :=0 ; 

TIME:=0 

Compute  E  =  the  sum  of  the  areas  of  the  circles,  i.e. 
n 

^   =  I      E(C^) 

i=l 

E(Cj^) 
STEP  1   Randomly  select  cycle  C.  with  probability  p.  =  

STEP  2_  Randomly  select  a  point  s  in  C. 

STEP  2  f : =0 

Repeat  until  a  C.  is  selected  such  that  seC^: 

BEGIN 


Randomly  select  C.  with  probability  1/n 
f :=f+l 


Test  if  s€  C. 

TIME  :=  TIME  +1 
END 
Let  U  :=  L'L 


STEP  4  number'of • trials :=number'of 'trials  +  1 

S  :=  S  +  U 

if  TIME  <  ^H   then  goto  STEP  1 
—        6^2  


STEP  5   Output  U*:^ 
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S 


number'of • trials 


End  of  algorithm  LMC . 

By  the  description  of  the  algorithm  it  is  clear  that  it   stops   in 
0[ J  steps. 

We  wish  to  prove  that  LMC  is  an  (e,5)  algorithm,  i.e.   that 

Prob{  |U*-U|  >  e  U|  <  5  i.e. 

Probj  IS  -  T-U|  >  eTU}  <  5  (INI) 

where  T  is  the  number* of • trials 


3.2  Notation  and  some  useful  remarks. 

Definition.   Let  t(T)  be  the  sum  f,  +  fj  +...+  f-n,  i.e.   the  sun  of  f's 
of  the  first  T  trials. 

Definition.   Let  U*(T)  be  the  value  of  U*   computed   if   the   algorithm 
makes  T  trials. 


Definition.    Let   U^^   be  the  estimate  of  the  area  of  the  union  for  the 
i"^^  trial. 


Definition.   Let  T(t)  be  the  number  of  trials  done  In  TIME  =  t 
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Deflnltton.   Let  U*(t)  be  the  output  of  algorithm  LMC  given  TIME  =  t. 

Note  that  the  details  of  STEP  2  are  as   follows:   The   probability 

that  the  point  s  will  fall  into  an  elementary  area  dA  of  C.  is  the  same 

for  all  elementary  areas  and  is  equal  to  dA  (R^^  =  radius  of  C^)  . 

TvR? 
Let  f^  be  a  random  variable  indicating  the  value  of  f  at   the   end 

of  the  i"^^  trial  in  LMC. 
Lemma  3 . 1 


(a)  The  mean  value  of  f j^  is  mean  (f^)    =  ^ 

(b)  The  mean  value  of  f?  is  mean  (f?)  <  i^lii, 


Proof 

(a)  Let  E(C.  ,s)  be  the  event  that  cycle  C.  is  selected  in  STEP  1 
of  the  i  trial  and  a  point  is  selected  in  C.  in  the  small  area  dA  (in 
C.  )  around  se  C-  . 

We  have 


E(C.)     , 
ProfJECC.  ,s)t  =  —J-   .  -L.  d 


.A 


Also 

Probjf  j^=k/E(C.  ,s)}  =   ( l-p(  success) )  )^~^  •  p(success) 


1.      /       N    cover(s)  „  ,  . 

where   p(success)  = .  We         conclude         that 


mean[fi/E(C   s)J  =     "  ,  , 
J       cover(s) 

But 


I        .';   ■ 
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E(C. )   , 
mean(fj^)  =  I  mean  (fl/E(C.  ,s))  •  ^ ^  dA 

alKCj  ,s)  ^   TtRj 


1  .e. 


mean(f ,)  =  —   J        L—.  dA 

IE-'    ,,    cover(s)  E(C.  ) 
over  all        \  /   v  ^  / 


circles  C. 


,,  ,    n   I  dA 

mean(f  j^)  =  -  J  ^r-^ 

'-          L  cover(s) 
over 


all  circles  C^ 

But   J         =  area  of  union  =  U.   So,  mean  (f.)  =  — •   QED 

cover(s)  '■  Z 

over 

all  circles 


(b)   mean  (f?/E(C^,s))  =  I     k'-C  l-p(success) )  ^    ^'pCsuccess) 

k=0 


2 
p  (success) 


(by  basic  properties  of  geometric  distributions) 


So 


^2 
mean(f|/E(C. ,s))  <  "" 


cover'^(s ) 
Hence 


9                  2n2     E(C.)      , 
mean(f2)  <   J        _JlL ^  .    L^  dA 


ill  (Cj  ,s)  cover^(s)         ''^^j 


trc  ■  rj;- '    f  idPi^ii 
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<:  2n^    r          dA     ,,           1^1, 
^  — — -   I       -. — ^  I  because  <  -. — r-l 

^   all  C.'s  ^°^^"^'^  cover2(s)    cover(s)^ 


So 


2 
mean  (f?)  <  1^—  QED 


Corollary  3.1 

Let  U.  be  the  random  variable  indicating  the  value  of  U  in  the  i 
round.   Then  mean  (U.)  =  U  and  mean  (U,)  <  2  E*U 


Proof 


U-  =  —  •  f  .  =»  mean  (U^ )  =  —  •  —  =  U. 
^   n     '■  ^    n    E 


Also 


"i  =  —  fi  ^  "^^^^   ^^0    *  -2    '  'T~ 
n  n 


(by  using  Lemma  3.1b) 

=>  mean  (U?)  <  2EU.  QED 


-15- 
3.3   The  proof  that  LMC  is  an  (e,5)  algorithm. 

Theorem  3.1   For  any  G,6e  (0,1) 

Prob{|U*(t)  -  U|  >  eU}  <  6    where 


en 
t=  ,  c  a  constant. 


Proof.   We  use  here  the  techniques  of  [Karp,  Luby  1983]. 
From  the  program  of  LMC 

f^  +  f2  +...+f'Y   =  t(T)    (T  =  total  number  of  trials)   i  .e , 
T'mean  (f^)  =  mean  (t(T)) 


c^n 

So,  if  t  =  for  some  c,  constant  then  we  expect  to  do 

6e^ 


^  ^  ^1  E  n  ^  '^l  E  ^  ^1" 


5c2  U  n   5^2  U    ^,2 


6e:2 


Z      *                      ^1  E 
number  of  trials  (because  _  <;  n    ))•   Let  k  be  equal  to .   Put  t' 

"  5g2  U 

Cj^n  o^ 

•  and  t"  =  t'(l+6)  where  B  is  a  constant,  to  be  determined. 
Let  us  estimate  the 


„   ,1  |U*(t")-U|  ^   , 
Probj >  ej 


This  is  equal  to 


E  <  n  •  max   E(C^)  <■   nU 
1 
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Prob{ lU*(t")-U|  >  eU  and  T(t")  <  k}  + 

Prob(|U*(t")  -  U|  >  eU  and  T(t")  >  k|  (EQA) 

But 

Probi|U*(t")  -  U|  >  eU  and  T(t")  <  k} 

<  Prob|T(t")  <  k\  (EQB) 

But   T(t")  <  k  *=*  fi  +  f2  +•••+  fk  >  ^"-      So 


ProblT(t")  <  k}  =  Probjfi  +...+  f^  >  t"; 


ProbI  ''^'"^'^   >  t'(l^B) 


So 


fl+...+f^    t' 


ProbJKt")  <  k\    <  Probll !l_ —-1    >  8} 

V  t   I  K.  ) 


,  25c^   , 

^ (by  Cheyshev  inequalitv). 


(This  is  because  k  =  — -  —  and  mean  ( )  =  :nean(f.)  =  —  and   t' 

c^n  5e2  U  k  ^  l 


b.JL 
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Hence  —  =  —  =  mean  )  I 

k    E  k      ^ 

From  (EQB)  we  finally  have 


2 
Prob||U*(t")  -  U|  >  eU  and  T(t")  <  k}  <  ?ll_  (EQC) 

CiB2 


Also 

Prob{|U*(t")  -  U|  J»  eU  and  T(t")  >  k} 


U,+. ..+U 
=  Prob{  3  r  >  k:  T(t")  >   k  and  \— i  -  U|  >  e  u} 


<  Prob{3  r  >  k:  1 L  -  U|  >  eU}  (EQD) 


Here   we   use    the    following   face   Kolmogorov ' s   inequality   (see 
[Feller, 57]  ) 

Let  Yj^,...,Yj^   be   independent   random  variables   with   the   same 

distribution  as  Y  and  let  mean(Y)  and  mean(Y^)  be  finite.   Let  a-(Y)  be 

the  variance  of  Y. 
I 
Let  Kj^   =      V   Y^  ,  then 
i=l 

Probf  Si:    1  <  Jl  <;  n  and  IR^  -  Jl'mean(Y)|  >  e«m«mean(Y)} 


<     o2(Y) 


2        7 
e  'ra  mean  (Y) 


By  applying  the  above  inequality  for  '^i  =  U^^  ,  i=l,...,T'  ,  m  =  T' ,   we 
get 


U^+. •.+U£ 


Probf  i' 


-  U 
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e-^T'U 


[To  verify  this,  put  mean  (U^^)  =  U  V^^-  and  Y  =  U^  (some  i)  and — 

•"^^"^"^  ^  <  2  -)  (T'  is  arbitrary).   Let  us  now  use  T'  =  2^«T  i=l,,../ 
Then  we  get 


U^+.  ..+Ujj 


P^  =  Probj  I:    \    <    i   <   N«2^  and 


-  U 


>  2e\ 


So  we  conclude 


Ui+. ..+U£ 


ProbI  Z:    N  2^"^  <    Z    <   N«2^  and 


-  U 


>  2e| 


2i:    1 


U  e2N2i 


(EQE) 


Now,  we  have 


U^+. •■+^1 


-  u 


Probj 


<  G  for  some  I    >   k\ 


n^+. ..+u^ 


<   5^   Prob{  i:    k«2'-  ^  <  £  <  k«2i  and 
i=l 


U       I    2' 


<   I   ^—^4   (by  EQE) 
i=l  "  iz/2)\   2^ 
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6e2  U 


-1  ^ 


—   (use  value  of  k) 


86 


That  is,  we  just  proved 


Probj  For  some  £  >  k. 


Ui+...+U( 


-  U 


>  eU}  < 


86 


(EQF) 


By  (EQA),  (EQB),  (EQC)  ,  (EQD)  ,  (EQF)  we  now  have 


CiS^   c^ 


Let  5'  =  —  [ +  8j.   Then,  if  c  =  ci(l+B),  we  have  proved   that   for 

any  e,  5'  S  (0,1) 


Prob{ |U*(t")  -  U|  >  G  U}  <  6' 


en 


where  t"  =  ,  c  a  constant  I.e.  LMC  is  an  (e,6')  algorithm.   QED 

6'e2 


■^-  i    t . 


'f ''  I    tp-   •  i   i..-  I  3:'i,e- 


.- <=-::}  X 
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4.   Applications  to  Intersection  problems 

One  can  use  the  linear  probabilistic  estimation  of  union  in  order 
to  check  (in  linear  time)  if  the  area  (volume)  of  intersection  of  n 
cycles  (spheres)  is  at  least  a«E  where  ae  (0,1)  and  E  is  the  sum  of 
their  areas  (volumes). 

Let  U  be  the  area  (volume)  of  the  union  of  the  n  objects.  Clearly 
U  =  E  means  the  objects  do  not  intersect  (we  are  leaving  out  zero  area 
intersections).  We  want  to  test  efficiently  if  E-U  >  aE  i .e .  if 
U  <  E(l-a).   We  run  the  union  estimation  algorithm  with  a  small  e  (e.g. 

e  =  — ,  k  >  10).   If  the  estimated  U*  <  E(l-a)  -  e  then  with  probability 

k 

at  least  1-5  the  area  of  intersection  is  at  least  a«E. 

The  above  technique  can  be  extended  to  objects  of  higher 
dimensions . 


V--J 


rJr 
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